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CHAPTER  I 
INTRODUCTION 


The  subject  of  this  thesis  is  the  development  of  a satellite 
communications  antenna.  The  fundamental  requirement  is  that  the 
radiation  pattern  of  the  antenna  be  broad  enough  to  cover  the 
entire  portion  of  the  earth  visible  to  it  in  the  optical  sense  as 
illustrated  in  Fig.  1-1.  The  desired  radiation  pattern  should  have 
a maximum  at  the  horizon  while  a local  minimum  is  permissible  mid- 
way between  the  horizon  look  angles.  This  minimum  is  allowable  due 
to  the  increased  path  distance  to  the  horizon  as  compared  to  the 
shortest  distance  between  the  satellite  and  the  earth's  surface. 

The  local  minimum  can  be  as  much  as  6 dB  below  the  maximum  radi- 
ation in  the  direction  of  the  horizon.  The  desired  bandwidth  for 
this  antenna  is  800-1000  MHz. 


Fig.  1-1.  Satellite  above  the  earth  with  a 60°  look  angle 
to  the  horizon. 
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The  UHF  satellite  antenna  under  investigation  is  shown  in  the 
photoqraph  in  Fig.  1-2.  The  antenna  consists  of  four  base-fed  _ 
monopoles  fed  in  phase  quadrature  which  produces  circular  pot  an  za~ 
“S  « the  «1s  if  the  antenna.  Three  major  factors  centre  ling  the 
radiation  pattern  produced  by  this  antenna  are  the  flare  angle  o 
the  blades,  the  spacing  between  the  blades,  and  the  location  of  the 
blades  on  the  groSnd  plane.  Results  for  various  fare  angles 
soacings,  and  locations  are  presented  in  Chapter  IV.  The  results 
in  Chapter  IV  are  calculated  using  the  Hybrid  Technique  which  com- 
bines the  Method  of  Moments  and  the  Geometrical  Theory  of  Dif- 
fraction [1,2]. 

The  Method  of  Moments  [3]  is  discussed  in  Chapter 
with  Galerkin's  Method  [3],  Point-matching  [3],  and  the  Reaction 
Method  [4,5]. 

An  introduction  to  the  Geometrical  Theory  of  Diffraction  [6] 
is  presented  in  Chapter  III.  Particular  attention  is  focused  on 
thePequations  for  the  diffracted  field  from  a perfectly  conducting 
wedae  The  *dae  diffracted  field  from  a ground  plane  which  is 
electrically  small  is  of  great  importance  to  the  radiation  pattern 

shape. 


CHAPTER  II 
METHOD  OF  MOMENTS 


A.  Introduction 

Associated  with  any  general  radiation  problem  there  is  an 
appropriate  integral  equation.  With  the  Method  of  Moments  it  is 
possible  to  approximate  the  solution  to  the  integral  equation  by 
a system  of  linear  algebraic  equations. 

Consider  a metallic  body  with  current  density  J_ on  it.  If  the 
body  is  a perfect  conductor,  then  the  boundary  conditions  require 
that  the  tangential  electric  field  be  zero  on  the  body's  surface 
as  follows: 


<2-1>  -tan  + -tan  = ^ 

$ 

itan  Is  a scattered  electric  field  generated  by  the  current  density 
-•  --Ian  1S  ^ie  tan9ential  component  of  an  incident  electric  field 
generated  by  a source  anywhere  on  or  outside  the  body.  Throughout 
the  remainder  of  this  chapter,  the  tan  subscript  will  be  dropped, 
and  it  will  be  understood  that  the  tangential  component  is  used 
unless  otherwise  specified. 

Equation  (2-1)  may  be  rewritten  as, 

(2-2)  -ES  = E1 


The  relationship  between  the  current  density  and  the  scattered 
field  may  be  written  as 

(2-3)  lop(J)  a - £S 

where  L0p  is  defined  as  an  operator  which  must  be  determined  for 
the  particular  problem  of  interest.  The  concept  of  linear  vector 
spaces  and  operators  can  be  used  to  write 

(2-4)  Lop(J)  = Ef 
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In  this  equation,  E.1  is  a known  excitation  function,  and  0_  is  an 
unknown  response  function  to  be  determined.  For  the  problem  considered 
in  this  thesis,  L0p  is  an  integral  operator  on  the  current  density  J. 
For  a given  source,  it  is  assumed  that  is  unique  and  may  be  expressed 
as 

(2-5) 

For  a given  problem,  the  space  of  functions  or  the  domain  of 
definition  of  the  operator  must  be  determined  along  with  the  functions 
resulting  from  the  operation  or  range.  If  J_  is  unique,  then  the 
operator  must  perform  a one-to-one  mapping  in  going  from  a subset 
containing  0_  to  one  containing  E1 . 

In  addition  to  finding  the  domain  and  range  of  the  operator, 
it  is  necessary  to  formualte  an  inner  product  or  reaction  [7], 

<J.,E>.  The  reaction  is  a scalar  quantity  which  satisfies  the 
following  conditions: 


(2-6)  <J,E>  = <E,J> 

(2-7)  <aJ^  + (5J.,E>  = + 

if 

(2-8)  <J*,0>  >0,  then  J t 0 

if 

(2-9)  <J*,J>  = 0,  then  J = 0 

where  a and  8 are  scalars  and  * denotes  the  complex  conjuagte.  The 
reaction  is  discussed  in  detail  in  Section  D. 
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There  are  four  basic  steps  in  solving  Eq.  (2-5)  by  the  technique 
known  as  the  Method  of  Moments: 

1.  is  expanded  in  a series  of  basis  functions , Jn 
which  span  J in  the  domain  of  L . 

2.  A set  of  testing  or  weighting  functions  are  defined 
and  a suitable  inner  product  determined. 

3.  The  inner  product  is  taken  and  a matrix  equation 
formed. 

4.  The  matrix  equation  is  solved  for  the  unknown,  the 
current. 

Once  the  current  is  found,  then  the  far-zone  field  pattern  and 
impedance  may  be  calculated  in  a straightforward  manner. 

Depending  upon  how  the  basis  functions,  testing  functions  and 
inner  product  are  chosen,  several  differnt  types  of  solutions  are 
obtained.  If  the  basis  and  testing  functions  are  chosen  to  be 
identical,  we  have  Galerkin's  method.  If  Dirac  Delta  functions 
are  used  for  testing  function,  we  have  a method  known  as  point- 
matching. If  we  use  the  Reaction  Integral  Equation  and  piecewise 
sinusoidal  basis  functions,  we  have  Richmond's  Reaction  method  [4], 
which  is  equivalently  a Galerkin  formulation.  The  remainder  of 
this  chapter  is  a discussion  of  each  of  these  methods.  The  Reaction 
Method  is  subsequently  used  to  obtain  the  results  in  Chapter  IV. 


B.  Galerkin’s  Method 

In  the  first  step  of  the  general  Method  of  Moments  solution,  the 
response  function  J in  Eq.  (2-4)  is  expanded  in  a series  of  basis 
functions  d]  >^2 * * on  a surface  S and  defined  in  the  domain  of 
the  operator  L^. 

(2-10)  J=  y I J 
— n n ~n 

The  In's  are  unknown  complex  coefficients.  The  solution  for  these 
gives  the  amplitude  and  phase  of  the  current  on  the  radiator.  Sub- 
stituting Eq.  (2-10)  into  Eq.  (2-4)  yields 

(2-">  Lop(n„4)  -i1 
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or,  using  the  linearity  of  L 


op 


> 


(M2)  I In  Lop  (4)  ' i!  • 

The  second  step  in  the  solution  is  to  define  a set  of  weighting 
or  testing  functions  111  >M2 » * * ' in  the  range  of  L . Next,  the  inner 
product  is  formed 


(2-13) 


1 1, 


<w  , 

-in’ 


^op  i^n> 


*9 


If  Wj,,  is  set  equal  to  then  the  formulation  of  the  problem  is 
known  as  Galerkin’s  method,  and  Eq.  (2-13)  becomes 


(2-14)  l I, 


•ftL, 

-fir 


L J > 
op  -fl 


= <J  , En> 
-m  — 


In  the  third  step,  the  inner  products  or  reactions  are  calculated 
and  the  matrix  equation  formed: 


(2-15) 


» Lnn  il > ^ » Lop  i-2> 


-1  * op  -1 
<~2’  Lop  —2" 


( 'A 


* ri 

<Ho » E > 


V v W 


In  compact  matrix  notation. 


(2-16)  [Z]  (I)  = (V) 


[Z]  represents  the  generalized  impedances,  (I)  the  generalized 
currents  and  (V)  the  generalized  voltages.  Solving  Eq.  (2-16)  for 
(I): 


(2-17)  (I)  = [zr1  (V) 

where  (I)  is  a column  matrix  containing  the  In's  in  Eq.  (2-10). 
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At  least  two  integrations  must  be  performed  for  each  impedance 
element  in  Eq.  (2-15).  One  of  these  integrations  can  be  avoided  by 
using  point-matching. 


C.  Point-Matching 

In  Galerkin's  method,  it  is  necessary  to  evaluate  the  generalized 
impedance  elements 


(2-18)  Zmn  = <J  , L (0  )> 
' mn  -m  op  -n 


If  numerical  integrations  are  involved,  they  may  be  tedious  and 
computationally  expensive.  This  problem  can  be  remedied  by  choosing 
Dirac  delta  functions  as  weighting  functions  in  Eq.  (2-15)  such  that 


^(S-S.,),  Lop(J1)<  <6(S-S-| ) , LQp(J2)>... 

(2-19) 

<<5  ( S-S.2 ) » Lq  p(^_i ) > 

• 

• 

• 

[J 

f<6{ S-^),  e]>\ 
= <6  (S_-S^)  ’ —2* 

W'V  eJ>  J 


where  is  the  distance  to  some  reference  point  and  S is  the  distance 
to  the  point  at  which  the  boundary  condition  is  bei reapplied.  Ei 
indicates  that  matching  is  being  done  at  point  1,  E£  indicates  that 
matching  is  being  done  at  point  2,  and  so  forth.  Any  further 
integrations  remaining  are  those  due  to  L . 
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The  use  of  Dirac  Delta  functions  is  physically  a relaxation  of 
the  electric  field  boundary  condition.  Nov/  the  boundary  condition 
is  only  applied  at  discrete  points  on  the  surface  and  not  continuously 
over  the  entire  surface,  thus  the  term  point-matching  describes  the 
process.  Accurate  solutions  are  obtained  if  a sufficient  number  of 
points  are  picked  to  match  boundary  conditions.  The  location  of 
these  points  is  also  important  for  accuracy.  Equispaced  points  yield 
good  results  for  scattering  problems  or  far-field  patterns  v/here  the 
observer  is  far-removed  from  the  radiator.  For  near- field  problems, 
however,  the  number  and  location  of  the  match  points  is  more  critical. 

D.  The  Reaction  Method 

Richmond  [4]  utilizes  the  reaction  integral  equation  (P.IE)  to 
derive  the  impedance  matrix  for  the  th in-wire  radiators  or 
scatterers  in  his  thin-wire  computer  program  [8].  The  unknown 
quantity  in  the  RIE  is  the  current  distribution  of  the  wire 
structure.  To  arrive  at  the  impedance  matrix  from  this  viewpoint,  it 
is  necessary  to  define  suitable  test  sources  and  expansion  modes. 

The  purpose  of  this  section  is  to  define  the  RIE  and  illustrate  its 
solution  for  the  impedance  matrix  with  the  use  of  test  sources  and 
expansion  modes. 

Consider  a closed  surface,  S,  of  a wire  structure  with  interior 
volumetric  region,  V.  There  is  an  external  source  ( J_-j , M-j ) in  the 
vicinity  of  the  wire  structure  which  generates  the  fields  (E_,  H ) . 

If  the  wire  structure  is  removed,  and  the  source  is  in  a homogeneous 
medium  (u,e  ),  then  the  incident  fields  (Ej , H-j)  are  generated.  The 
scattered  fields  are  defined  as  the  difference  between  the  fields 
(£,  H)  and  (E^ , H. ) with  the  wire  structure  present. 

(2-20)  ^ = E - E. 

(2-21)  Hg  = ii  - H i 

These  are  time-harmonic  fields  at  the  same  frequency  with  the  time- 
dependence  term  ejut  suppressed. 

From  the  surface-equivalence  theorem  of  Schelkunoff  [9],  the 
following  surface-current  densities  can  be  assumed  on  S: 


(2-22)  Jg  = n x H 
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(2-23)  M = E x n 
-s  — 


The  fields  (JE,  Hj  interior  to  S can  be  assumed  to  be  zero,  while  the 
exterior  fields  (£.  H)  remain  unchanged,  allowing  the  wire  structure 
to  be  replaced  by  homogeneous  medium  (u,e).  The  unit  vector  n is 
directed  outward  from  the  surface. 

J and  Ms  generate  the  fields  (Eg,  Hs)  exterior  to  the  area  that 
once  contained  the  wire  structure.  Since  E_  and  H are  zero  interior 
to  the  region  once  containing  the  wire  structure,  from  Eqs . (2-20) 
and  (2-21)  these  interior  fields  are  -£.  and  -H . . 

The  Carson  reciprocity  theorem  [10]  is  defined  as 


(2-24) 


(J2  • E,  - H2  • H1)dv2 


(ill  > Mi)  ar*d  (M2»  M2)  are  source,  or  impressed,  current  densities  and 
(Ml’  Ml)  are  fields  from  source  1 radiating  in  the  presence  cf 
source  2,  and  vice  versa  for  (jip,  H?) . Either  side  of  Eq.  (2-^4)  is 
known  as  a reaction:  c 


(2-25)  <1,  2> 


• E2  - M,  • HjJdV, 


If  reciprocity  conditions  are  satisfied,  then 


(2-26)  <1,  2>  = <2,  1> 


Equation  (2-26)  can  be  expanded  +0  yield  Eq.  (2-24)  (the  Carson 
reciprocity  theorem).  Reaction  is  considered  as  a measure  of  the 
coupling  between  two  sources,  1 and  2. 

Going  back  to  our  original  problem,  a test  source  (Jjn,  Mm)  is 
placed  in  the  interior  region.  The  reaction  of  this  test  source 
with  the  other  sources  is  desired.  Applying  the  Carson  reciprocity 
theorem  and  integrating  over  the  surface  of  the  test  source: 


(2-27)  W.  • E,  - V H,)*  - - 


( V " V -i'ds 
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where  (£$,  Jis)  are  the  fields  generated  by  (Js,  M^),  and  the  origin 
of  ( Ej  ,~hi ) was  discussed  previously.  Eq.  (2-27Tis  one  form  of 
the  reaction  integral  equation  ( RI E) . Eq.  (2-27)  states  that  the 
interior  test  source  has  zero  reaction  with  the  other  sources. 
Another  form  of  the  RIE  may  be  obtained  from  Eq.  (2-27)  and  the 
reciprocity  theorem: 


(2-28) 


“s 


• Hm)ds  + 


f[[ 

JJ] 


(J.  • Em  - M.  • Hm)dv  = 0. 


(Em,  Hm)  are  the  fields  generated  by  the  test  source.  The  volume 
integral  is  over  that  region  which  used  to  be  occupied  by  the  wire 
structure.  Electric  test  sources  are  used  with  Eq.  (2-2C)  to  determine 
the  surface  current  distributions  83  and  Mr  for  thin-wire  problems. 

With  electric  test  sources  applied  to  Eq.  (2-28),  a less  general  form 
of  the  RIE  is  found  known  as  the  electric-field  integral  equation 
(EFIE). 


In  the  wire  structure,  let  each  segment  have  a circular  cylindrical 
surface.  A circular-cyl indrical  coordinate  system  with  unit  vectors 
(p,  $,  z)  is  set  up.  The  integral  equation  is  simplified  by  assuming 
the  wire  radius  a to  be  much  less  than  the  wavelength  A,  and  that  the 
length  of  the  wire  is  much  greater  than  its  radius.  Therefore,  if 
the  wire  is  thin,  the  integrations  over  the  ends  of  the  wire  can  be 
neglected  along  with  the  circumferential  component  of  the  surface- 
current  density  0,.  Furthermore,  the  axial  component  of  the  current 
density  Jz  can  be^assumed  independent  of  <j>.  By  these  assumptions, 
the  current  density  on  the  thin  wire  structure  is  related  to  the 
total  current  by 


(2-29) 


j tz)  = z Hz).  _ 

Vz'  2 it  a 


K z) 

2tt  a 


z is  a metric  coordinate  measuring  position  along  the  wire  axis  and 
I(z)  is  the  total  current  including  conduction  and  displacement.  If 
one  wire  is  within  a few  diameters  of  another  or  the  wire  is  bent 
to  form  a small  acute  angle,  then  a more  elaborate  formulation  than 
Eq.  (2-29)  is  required. 

If  the  wire  structure  is  perfectly  conducting,  then  Ms  is  zero. 
Substituting  Eq.  (2-29)  into  the  RIE  (Eq.  (2-28)),  Eq.  (2-28)  reduces 
to 


(2-30) 


•L 

0 


Kz)  En2'  dz  = V 


m 
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where  L is  the  entire  length  of  the  wire. 


(2-31) 


• Em  - M.  • Hm)  dv 


and 

(2-32)  = 2~~~  z • Em  d* 


The  sinusoidal  reaction  formulation  for  thin  wires  is  based  on 
Eq.  (2-30).  The  known  quantities  in  this  equation  are  £_m  and  V™. 

The  unknown  is  the  current  distribution  I(z).  To  solve  for  I(z), 
suitable  test  sources  and  expansion  modes  need  tc  be  defined. 

A good  choice  for  a test  source  for  a perfectly  conducting 
thin  wire  is  a filamentary  electric  dipole  with  a sinusoidal  current 
distribution.  The  piecewise-sinusoidal  function  is  similar  to  the 
natural  current  distribution  of  the  perfectly  conducting  thin  wire. 
Also,  the  sinusoidal  dipole  represents  a finite  electric  line  source 
in  a homogeneous  medium  with  simple  closed-form  expressions  for  the 
near-zone  fields,  while  the  mutual  impedance  between  two  such  dipoles 
can  be  represented  in  terms  of  exponential  integrals. 

A V dipole,  with  different  arm  lengths  and  terminals  at  the 
vertex,  makes  a typical  test  source.  The  current  is  a maximum  at  the 
terminals  and  decreases  sinusoidally  to  zero  at  the  endpoints.  The 
terminal  current  is  one  ampere,  and  the  current  distribution  has  a 
slope  discontinuity  at  the  terminals.  The  linear  test  dipole  along 
with  its  current  distribution  is  shown  in  Fig.  2-1. 


Fig,  2-1.  A linear  test  dipole  and  its  sinusoidal  current  distribution 
The  endpoints  are  at  z-|  and  Z3  with  terminals  at  Z2- 
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The  current  distribution  J.(z)  = £.(z)  can  be  represented  as 


(2-33) 


z P,  sin  k(z-z,)  z P9  sin  k(z--z) 
Hz)  = — U— „ l + L 3 


sin  kiz^'Z-] ) 


sin  kflTpip 


where  P]  = 1 on  the  interval  Z]<z<Z2  and  zero  elsewhere  and  P2  = 1 
for  Z2<z<Z3  and  is  zero  elsewhere,  k is  the  complex  propagation  con- 
stant of  tne  homogeneous  exterior  medium: 


(2-34)  k = jW  ijc 


This  value  of  k is  necessary  for  the  assumptions  made  earlier 
regarding  the  test  source. 

To  simplify  the  integration  in  Eq.  (2-32),  the  test  dipole  is 
placed  on  the  wire  axis  rather  than  elsewhere  in  the  interior  region 
of  the  wire  structure. 

A practical  problem  may  require  many  test  dipoles  located  at 
different  positions  along  the  wire  axis  to  form  an  overlapping  array. 

A particular  test  dipole  nt  will  generate  the  fields  (JP,  Jf11)  in  the 
homogeneous  medium.  Each  of  the  test  dipoles  must  satisfy  Eq.  (2-30). 
Also,  each  dipole  must  radiate  at  the  same  frequency  as  the  wire 
structure's  true  source.  Thus,  enforcing  Eq.  (2-30)  requires  each 
test  dipole  in  the  array  to  have  the  correct  reaction  with  the  true 
source.  If  there  are  N test  dipoles,  then  Eq.  (2-30)  represents  a 
system  of  N simultaneous  integral  equations. 

The  current  distribution  on  the  wire  structure  may  be  expanded 
in  a finite  series: 

N 

(2-35)  I(z)  = l In  Fjz) 

Fn(z)  is  a normalized  expansion  function  equal  to  the  test  dipole 
current  distribution  in  Eq.  (2-33).  Since  a single  expansion  function 
relates  to  only  a two-segment  portion  rather  than  the  entire  wire 
structure,  these  functions  are  commonly  referred  to  as  subsectional 
bases.  The  In's  in  Eq.  (2-35)  are  complex  constants  which  represent 
samples  of  the  current  function  I(z).  The  sinusoidal  bases  will 
resemble  the  triangular  bases  of  a piecewise-1 inear  model  as  il- 
lustrated in  Fig.  2-2,  if  the  wire  segments  are  electrically  short 
in  reference  to  wavelength.  Also  illustrated  in  Fig.  2-2  is  the 
current  distribution  I(z)  and  its  two-mode  approximation  1 1 ( z ) . 
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F,(z)  F 2 (2.) 


Fig.  2-2.  The  expansion  functions  F-j (z)  and  F2(z),  the  current 
distribution  I(z)  and  the  two-mode  approximation 
I'(z). 


Notice  the  slope  discontinuities  of  the  piecewise-sinusoidal 
expansion.  These  discontinuities  occur  at  the  generators,  lumped 
loads  and  wire  corners.  For  all  normal  calculations,  the  calculated 
samples  In  will  be  accurate,  and  the  corresponding  piecewise- 
sinusoidal  current  distribution  I * (z)  will  be  satisfactory. 

By  substituting  Eq.  (2-35)  into  Eq.  (2-30),  the  following 
set  of  simultaneous  linear  algebraic  equations  is  obtained: 

N 

(2-36)  l I Z = V where  m = 1 ,2 , • . -N 

n=l 


and  where 


(2-37) 


inn 


Fn(z)  ^ dZ 
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The  integral  in  Eq.  (2-37)  has  limits  which  extend  over  the  two 
segments  in  the  range  of  expansion  function  Fn.  For  ease  of  computer 
manipulation,  Eq.  (2-36)  may  be  expressed  in  matrix  notation: 


<2-38>  tZJ  <‘n>  * <V„,> 

where  [Znm]  is  a square  matrix  and  (In)  and  (Vm)  are  the  column 
current  and  voltage  matricies,  respectively. 

[Zmn]  is  the  open-circuit  impedance  matrix  for  the  wire-structure 
where  the  first  subscript  m indicates  the  row  and  the  second,  n, 
the  column.  Zmn  also  represents  the  mutual  impedance  between  test 
dipole  m and  expansion  mode  n. 
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CHAPTER  III 

THE  GEOMETRICAL  THEORY  OF  DIFFRACTION 


A.  Introduction 

The  Method  of  Moments  is  usually  referred  to  as  a low  frequency 
technique  as  its  use  is  usually  applied  to  bodies  small  in  terms 
of  wavelength.  The  geometrical  theory  of  diffraction  (GTD),  however, 
may  be  applied  to  bodies  which  are  arbitrarily  large  in  the 
electrical  sense  and  is  thus  referred  to  as  a high  frequency  technique. 
GTD  employs  rays  in  a systematic  way  to  describe  the  field  scat  'red 
by  specific  parts  of  a body,  e.g.,  edges,  tips,  or  corners.  Only 
perfectly  conducting  bodies  in  an  isotropic,  homogeneous  medium  will 
be  considered  in  the  discussion  that  follows. 

The  basic  problem  in  GTD  is  the  two-dimensional  case  of  the 
diffraction  of  an  electromagnetic  wave  from  the  edge  of  a perfectly 
conducting  wedge.  Since  this  edge  diffracted  wave  behaves  like  a 
cylindrical  wave  radiating  from  the  edge,  the  contribution  from  an 
edge  element  on  a three-dimensional  body  may  be  approximated  by 
assuming  the  element  is  on  the  edge  of  a wedge  extending  to  infinity. 
Wedge  diffraction  can,  therefore,  be  applied  to  three-dimensional 
geometries  having  finite  edges. 

GTD  was  originally  based  on  plane  wave  diffraction  coefficients 
when  applied  to  wedge  diffraction.  However,  a different  formulation 
for  the  wedge  diffraction  is  substituted  depending  upon  whether  the 
incident  wave  is  cyl indirical , conical,  or  spherical.  Rudduck  [6] 
showed  that  in  the  treatment  of  antennas,  it  was  necessary  to  utilize 
the  diffraction  of  cylindrical  waves. 

B.  Wedge  Diffraction 

Consider  the  problem  of  a plane  wave  incident  upon  a wedge 
with  angle  (WA)  equal  to  (2-n)n  as  shown  in  Fig.  3-1.  At  the 
observation  point  P(r,$)  the  z-component  of  the  total  field  is 
represented  by  u(r>(j>) 

(3-1)  u(r,$)  = V(r,  <j>  - $'}  ± V(r,  $ + <{.'). 


The  plus  sign  is  used  if  the  electric  field  is  polarized  perpendicular 
to  the  perfectly  conducting  edge,  and  the  minus  sign  is  used  if  the 
electric  field  is  parallel  to  the  edge.  V(r,  «{>-(+>')  represents  the 
incident  field  and  V(r,  $ + $')  represents  the  reflected  field. 
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Fig.  3-1.  Conducting  wedge  showing  plane  wave  incidence. 


These  field  quantities  may  be  divided  into  the  geometrical  field 
V*  and  the  diffracted  field  V^, 


(3-2)  V(r , b)  = V*(r,  e)  + Vg(r,  B). 
The  geometrical  field  is  given  by 


(3-3) 


V*(r,  e)  = 


eJkr  cos  (b  + 2u  n N) 

! 

0 


if  - u<8+  2tt  n N < 7i 
otherwise 


for  N = 0,  ±1 , ±2, • • • . 
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For  (3  = <ji  - 4)',  Eq.  (3-2)  gives  the  incident  field.  The  incident 
field  will  be  seen  by  any  observer  in  regions  I and  II  of  Fig.  3-1. 

This  region  is  referred  to  as  the  illumination  region  (0<4><tt+4>  1 ) 
for  the  incident  field.  At  8 = n+  the  incident  shadow  boundary 
(ISB)  is  encountered.  Beyond  this  angle  is  the  shadow  region 
( n+<j> 1 <4  <itn).  Any  observer  in  the  shadow  region  will  not  see  an 
incident  field.  For  {$  = <j>  + <j>' . Eq.  (3-2)  yields  the  reflected 
field.  The  illumination  region  for  the  reflected  field  is  region  I 
of  Fig.  3-1  (0<4><ir-t{> 1 ) . At  6 = u - <j>'  is  the  reflection  shadow 
boundary  (RSB).  Any  observer  beyond  the  RSB  and  outside  of  region  I 
is  in  the  shadow  region  for  the  reflected  field  and  will  not  see 
that  component  of  the  field. 

As  indicated  by  Eq.  (3-2),  each  geometrical  field  has  a dif- 
fracted field  Vg( r , g)  associated  with  it.  The  diffracted  field 
takes  into  account  the  diffraction  phenomenon  at  the  edge  of  the 
wedge  and  combines  with  the  geometrical  field  to  eliminate  dis- 
continuities on  the  shadow  boundaries. 

The  diffraction  problem  from  a conducting  wedge  was  first  solved 
by  Sommerfeld  [11].  His  solution  was  in  the  form  of  a contour  integra- 
tion which  proved  inconvenient  for  computational  purposes.  He  later 
derived  a more  practical  solution,  but  it  was  only  valid  for  large 
values  of  kr(l  + cos  g) . The  form  of  this  solution  is 


(3-4) 


VB(r,  6)  = - 


-j(kr  + fj 
J 2*  kr 


cos  £ - cos  f- 
n n 


Pauli  [12]  developed  a practical  formula  for  the  field  diffracted 
by  a conducting  wedge  which  turned  out  to  be  more  accurate  than 
Sommerfeld 's  solution,  especially  in  the  vicinity  of  the  shadow 
boundaries , 


(3-5) 


. His  solution  i 

is : 

Vp(r,  g) 

i- 

_ 2 e 4 

(l  sin 

o' 

V" 

eJkr  cos 

8 

•CD 

jTfcn 

cos  6) 

cos 


i 


cos  £ - cos  & 
n n 


-Jt 


dT  + [Higher  Order  Terms] 


The  higher  order  terms  may  be  neglected  for  large  values  of  kr. 
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A diffraction  coefficient  D(p)  designated  by  Keller  [13]  may  be 
related  to  the  function  by 

j kr 

(3-6)  VB(r,  6)  = DU)  ° 


In  1967  a formulation  introduced  by  Hutchins  and  Kouyoumjian 
[14]  yielded  far  better  accuracy,  particularly  near  the  shadow  bound- 
aries and  for  r<x  than  that  obtained  from  the  Pauli  formulation 
(Eq.  (3-5)).  Hutchins'  and  Kouyoumjian's  formulation  is  given  by, 


(3-7)  Vb(L,  6,  n)  = I_w  (L,  p,  n)  + I+w  (L,  p,  n), 


where 

(3-8) 


where 


-j(kL+^-) 

(Ll  6>  n)  = F cot(^ 


JkLa 


4 2 


jTtE 


e Jt  dt  + [Higher  Order  Terms] 


(3-9)  a = 1 + cos  (p  - 2 nnN). 


N is  a positive  or  negative  integer  or  zero,  whichever  best  satisfies 
the  aquations, 


(3-10)  2 nnN  - 6 = - n , for  I 

and 

(3-11)  2 nnN  - 6 = +it  , for  I 


L is  a distance  parameter  to  be  defined  later. 
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Consider  a source  field  E^(s')  located  at  s ' ( p 1 ,<}>')  as  depicted 
in  Fig.  3-2.  The  source  may  be  either  electric  or  magnetic  and  may 
generate  an  incident  wave  on  the  edge  of  the  wedge  which  is  either 
plane,  cylindrical,  conical,  or  spherical.  The  diffracted  vector 
field  at  the  point  s Cp » 4>)  may  be  represented  in  terms  of  a dyadic 
diffraction  coefficient.  The  diffracted  fields  may  be  represented 
in  a more  compact  form  if  they  are  expressed  in  terms  of  a ray-fixed 
coordinate  system  centered  at  the  point  or  points  of  diffraction  QF 
as  suggested  by  Kouyoumjian  and  Pathak  [15].  t 

The  diffracted  field  may  be  written  as 
(3-12)  Ed  (s)  = Ei(QE)  • De  (s,  I)  A(s)e'jks. 


The  relationships  between  the  orthogonal  unit  vectors  is  shown  in 
Fig.  3-2.  Writing  the  diffracted  field  in  a more  convenient  form 
involving  the  VB  function, 


(3-13) 


<!<*> 

-VB 

0 

EX> 

_ Ef(s)  _ 

% 

0 

-v+ 

B J 

. EX>. 

JT  ejkL 
Sln 


A(s)e“jks 


where 


(3-14) 


VB(L,g",  n)  + Vb(L,  S+,  n). 


The  minus  sign  is  used  for  an  electric  field  polarized  parallel  to 
the  edge  (soft  boundary  condition),  and  the  positive  sign  is  used 
for  the  electric  field  polarized  perpendicular  to  the  edge  (hard 
boundary  condition).  g+  = <}>  +<$>'  where  the  minus  sign  is  used  for 
the  incident  field  and  the  positive  for  the  reflected  field.  The 
ray  divergence  factor  or  spatial  attenuation  factor  A(s)  is  defined 
as 


(3-15) 


A(S) 


1 

F 


= < 


s(s^s) 


for  plane,  cylindrical  and 
conical  wave  incidence 


Xr  spherical  wave  incidence 


The  distance  L 


(3-16)  L = 


parameter  is  defined  as 


s sin2 


P-..P- 
p 


+ p' 


2 

s's  sin 


s + s' 


for  plane  wave  incidence 
for  cylindrical  wave  incidence 

for  conical  and  spherical  wave 
incidence 


Through  the  introduction  of  L by  Pathak  and  Kouyoumjian,  it  is 
possible  to  use  the  diffraction  formulation  of  Eq.  (3-13)  for  near- 
zone cases.  Experience  has  shown  that  accurate  results  may  be 
obtained  for  kL  > 1 . This  formulation  will  prove  useful  in 
Chapter  IV  when  we  consider  a source  near  a conducting  wedge. 


CHAPTER  IV 


A HYBRID  TECHNIQUE  COMBINING  THE  METHOD  OF  MOMENTS 
AND  THE  GEOMETRICAL  THEORY  OF  DIFFRACTION 


A.  Introduction 

Consider  the  problem  of  a monopole  near  the  edge  of  a perfectly 
conducting  wedge  as  shown  in  Fig.  4-la.  To  determine  the  current 
distribution  on  the  monopole,  it  is  necessary  to  change  Eq.  (2-38) 
to  include  the  effect  of  the  edge.  This  modified  equation  for  the 
impedance  matrix  is  given  as  (1,2] 


<4-’>  &L*  <‘n>  ’ <vm> 

where  [Z^n]  is  the  new  impedance  matrix.  (I1)  is  the  modified 
column  current  matrix  due  to  the  effect  of  the  edge.  An  element 
of  [Z^  ] may  be  expressed  as 


(4-2)  Z'  = l + aZ 
' mri  mn  mn 


Z.im  is  the  original  impedance  maf  ix  from  Eq.  (2-37)  and  aZ  is 
tne  additional  term  due  to  the  effect  of  the  edge.  !l ' 

To  gain  some  physical  insight  into  the  problem,  consider 
Fig.  4-lb.  The  monopole  is  segmented,  and  on  each  segment  the 
current  is  assumed  to  be  constant.  Zmn  is  interpreted  physically 
to  be  due  to  the  z component  of  the  electric  field  at  segment  m 
from  the  unit  current  passing  through  segment  n.  AZ(i1n  is  defined 
as  the  reaction  between  a field  diffracted  from  the  edge  of  the 
wedge  due  to  a unit  current  flowing  in  segment  n„  E<j,  »,  and  the 
current  density  in  segment  m J as  follows:  ^ ' 

(4-3)  AZ  = <0  , Ed,  \> 
v 1 inn  m*  n(u 

The  geometry  of  the  situation  is  shown  in  Fig.  4-lc.  Ed/X\  is  found 
by  applying  geometrical  theory  of  diffraction  techniques.  To  calcu- 
late AZ,nn,  the  near-field  incident  upon  the  edge  due  to  segment  n is 
found  using  the  method  of  moments.  The  component  of  this  field 
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Fig.  4-la.  Monopole  on  a conducting  wedge. 


r 


Fig.  4-lb.  Monopole  on  a conducting  wedge  with  segment  m 
receiving  radiation  from  source  segment  n. 


Fig.  4-1 c.  Monopole  on  a conducting  wedge  with  segment  m 
receiving  diffracted  radiation  from  source 
segment  n. 


which  is  perpendicular  to  the  edge  and  to  the  direction  of  propaga- 
tion is  the  term  E^(Qe)  in  Eq.  (3-13).  This  component  is  found  from 
the  total  field  by: 


(4-4)  e]  (Qe)  = E^  sin  + e’  cos  ' 

E^1  and  E^  are  the  p and  z circular  cylindrical  coordinate  components 
of  the  electric  near  field  due  to  the  current  on  segment  n of  the 
monopole.  The  field  diffracted  from  the  edge  is  now  calculated  by 
applying  Eq.  (3-13) : 


(4-5) 


(s)  = - 


+ * 
V E' 
VB  -t 


(Qc) 


JT  ejkL 
sin 


A(s)  e 


-jks 


S0  = 90°,  therefore,  sin  6o  = 1 • Eor  spherical  wave  incidence 
( s.  \l/2  s.  s 

A(s)  "L(s,  and  L = ^r-  . Equation  (4-5)  may  now  be 

written  as 
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(4-6)  Ej(s)  = -Vg  e}  (Qe) 


s'  s 

1/2 

s 1 

s'  + s 

s(s  + s') 

1/2 


jk 


— 1 
+ sj 


,-Jks 


Simplifying, 
-d 


(4-7)  Ei  (s)  = -V*  El  (Qe) 


jk 

e l 


ss 


s + s 


r -s 


Equation  (4-7)  is  then  substituted  for  E*j/,\  in  Eq.  (4-3),  and  the 
reaction  integral  is  carried  out: 


(4-8) 


m 


(s)  cos  $ d£m 


where  the  integral  is  over  t!ie  length  of  segment  m. 

For  each  of  the  impedance  elements,  Zmr,,  of  Eq.  (2-38),  the 
modified  impedance  element  Z^n  is  determined  by  summing  Zmn  with  the 
corresponding  AZmn.  Z^n  is  then  substituted  for  Zmn  in  Eq.  (2-38) 
which,  when  solved,  will  result  in  a new  current  (I1). 

B.  Far-Zone  Field  Pattern  Calculation 

The  far-zone  field  patterns  are  determined  by  using  the  hybrid 
technique  and  a two-point  diffraction  method  as  illustrated  in  Fig. 
4-2a  to  approximate  the  true  physical  situation  depicted  in  Fig.  4-2b. 
The  two-point  diffraction  method  assumes  that  the  ground  plane  extends 
to  infinity  in  the  x direction,  allowing  us  to  ignore  the  contribution 
of  the  diffracted  fields  from  those  edges  and  the  corners.  The 
true  physical  situation  will  have  diffracted  field  contributions 
from  all  four  sides  and  corners,  but  the  contributions  from  points  1 
and  2 in  the  o-plane  will  dominate.  The  sum  of  X,  and  Xp  in  Figs. 

4-2a  and  b equals  28  inches. 

The  patterns  are  calculated  by  summing  the  direct  electric  field 
due  to  the  new  current  (I1)  on  segment  n with  the  incident  near  field 
from  segment  n diffracted  from  points  1 and  2 into  the  far  zone, 

(4-9)  E1 F = Er 1 + E J + Ep 
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sail kjjfe*. 


where  E_n  is  the  direct  field  from  segment  n,  and  E-?  are  the 
fields  diffracted  from  points  1 and  2,  respectively,  and  Ej  is  the 
total  far-zone  field  due  to  segment  n.  The  geometry  is  depicted 
in  Fig.  4-2c.  At  a given  angle  of  o,  the  total  fields  for  each 
of  the  monopole's  segments  are  summed  to  give  the  particular 
far-zone  field  at  that  angle. 


Fig.  4-2a.  Two-point  diffraction  approximation  for  determining 
far-field  patterns  for  a single  monopole  on  a ground 
plane.  X,  + Xo  = 28". 


180° 

z 


Fig.  4-2b . True  physical  situation. 


X1  + X2 


The  direct  electric  field  from  segment  n is  calculated  via  the 
Method  of  Moments,  while  the  diffracted  fields  from  the  edges  are 
found  using  the  Geometrical  Theory  of  Diffraction.  The  same  procedure 
applied  before  in  this  chapter  to  find  the  modified  impedance  matrix 
is  used  here.  However,  instead  of  finding  the  edge  diffracted  field 
at  another  segment  m on  che  mor.oDole  due  to  a segment  n,  the  edge 
diffracted  field  in  the  far  zone  is  desired.  Applying  Eq,  (3-13), 
we  find  the  perpendicular  component  of  the  diffracted  field  to  be, 


(4-10)  eJ  (s)  = -V*  Z[  (Qf) 


s' 

s + s' 


s'  s 

s + s'J  -jk« 


Jkx  sin  o 


W 


Fig.  4-2c,  Contributions  to  the  far-zone  field  include 
the  direct  field  E_n  and  the  diffracted  fields 
E?  and  eR. 


-1 


—2” 


where  x is  the  distance  from  the  edge  to  the  point  where  the  phase 
is  referenced,  and  the  term  ejkx  sino  puts  the  diffracted  field  in 
the  proper  phase  with  the  direct  field.  In  the  far  zone  s is  much 
much  greater  than  s';  therefore,  Eq.  (4-10)  becomes 


(4-11) 


El  (s)  = -V*  E}  (Qe)  s1  e 


Jks’ 


ejkx  sino  e 


jks 

s 


e"jks 
The  term  — - — 


is  suppressed  in  the  far  zone  reducing  Eq.  (4-11)  to 


(4-12) 


El  (s)  = -Vg  E|  (C^)  s' 


eJks 1 ejkx  sin  o 


Since  the  electric  field  is  polarized  perpendicular  to  the 
edge  of  the  ground  plane  Vjt  is  used.  The  significance  of  the 
+ and  - signs  was  discussed  in  Chapter  III. 
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Fig.  4-2d.  Geometry  for  the  second  order  diffracted  fields. 


To  take  care  of  discontinuities  occurring  in  the  patterns  at 
the  edges  of  the  ground  plane,  it  was  necessary  to  include  the 
second  order  diffracted  fields  from  edges  1 and  2.  The  second  order 
diffracted  field  is  an  edge  diffracted  field  whose  incident  field 
is  the  diffracted  field  from  the  opposite  edge  due  to  a particular 
monopole  segment  n.  The  geometry  is  shown  in  Fig.  4-2d.  There  are 
two  steps  in  finding  the  contribution  of  the  second  order  diffracted 
fields  in  the  far  zone.  First,  the  near-field  diffraction  from 
one  edge  and  incident  upon  the  opposite  is  calculated.  For  this,  a 
form  like  Eq.  (4-7)  is  used, 


(4-13)  E“  (s)  = -V*  e[  (Qe) 


jk 


s s' 
s + s' 


-s 


s + s' 


Two  fields  must  be  calculated  using  Eg.  (4-13).  One  field  travels 
across  the  top  of  the  ground  plane  (E^op)  and  the  other  travels  along 
the  underside  (Efower)  as  shown  in  Fig.  4-2d.  Since  each  of  these 
fields  are  grazingincident  fields  (direction  of  propagation  parallel 
with  ground  plane),  they  are  each. divided  by  a factor  of  two  and 
substituted  into  Eq.  (4-12)  for  EMQfr).  The  two  results  are  summed 
to  yield  the  second  order  diffracted  field  in  the  far  zone.  The 
second  order  diffracted  fields  are  calculated  for  both  edges 
(E?D  and  E|u)  and  then  summed  in  Eq.  (4-9)  to  yield  a new  expression 
for  the  total  far- zone  field  as  follows, 
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(4-14) 


E T = E n + E f + E ° + Ef  + £2° 


With  these  expressions,  the  far-zone  field  patterns  of  one  monopole 
or  two  monopoles  on  a ground  plane  are  accurately  approximated. 


C.  Single  Monopole  on  a Ground  Plane 

The  far-field  pattern  of  a quarter  wavelength  monopole  calculated 
using  the  two-point  diffraction  approximation  is  shown  in  Fig.  4-3a. 
The  X]  and  X2  values  in  Fig.  4-2a  are  each  three  wavelengths  as  in- 
dicated by  the  computer  drawn  sketch  in  Fig.  4-3a.  This  pattern  is 
compared  with  the  measured  pattern  of  the  actual  situation  and  the 
calculated  pattern  of  Lopez  [16]  in  Fig.  4-3b.  Lopez's  patterns  are 
those  of  a stub  on  a circular  ground  plane  with  a three-wavelength 
radius.  Lopez  also  used  a two-point  diffraction  approximation. 

His  two-point  diffraction  method  is  valid  everywhere  except  near 
the  caustic  which  exists  in  his  problem.  By  studying  the  measured 
and  calculated  patterns  in  Fig.  4-3b,  it  is  evident  that  this  is 
indeed  the  case.  From  0 to  about  10  degrees  in  Fig.  4-3b,  the  shapes 
of  the  calculated  and  measured  patterns  differ.  It  is  seen  from  Figs. 
4-3a  and  b that  the  two-point  diffraction  method  is  a reasonable 
approximation  for  the  true  physical  situation  as  the  two  points  are 
tlie  dominate  diffraction  mechanisms,  and  the  only  other  contributor 
to  the  far-field  pattern  is  the  monopole. 

In  Figs.  4-4a  through  4-6a  the  far-field  patterns  are  calculated 
using  the  hybrid  technique  and  two-point  diffraction  at  14  inches  to 
either  side  of  the  1/4a  monopole.  The  a parts  of  Figs.  4-4  through 
4-6  are  to  be  compared  with  the  b parts  which  are  far-field  patterns 
of  a stub  on  a four-sided  ground  plane  calculated  by  Burnside's  GTD 
computer  program  [17].  The  reason  for  making  the  comparison  with 
the  pure  GTD  program  is  to  provide  an  independent  check  on  the 
validity  of  our  computer  program. 

Burnside's  program  uses  infinitesimally  short  current  sources 
which  have  amplitudes  and  locations  so  as  to  simulate  a quarter 
wavelength  monopole  with  a sinusoidal  current  distribution.  His 
program  also  takes  into  account  the  diffraction  from  the  four 
corners  and  all  four  sides.  However,  in  this  comparison,  the  two 
opposite  sides  are  set  at  distances  much  greater  than  the  other 
two  so  as  to  approximate  Fig.  4-2a  and  so  that  the  points  on  the 
two  closer  sides  will  dominate  the  diffracted  field  contribution. 
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Fig.  4-3b.  Radiation  pattern  of  a stub  on  a circular  ground 
plane. 
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4-4b.  Calculated  by  Burnside's  method. 
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Fig.  4-5b,  Calculated  by  Burnside's  method 


Fig.  4-6a, 


~ig.  4-6b.  Calculated  by  Burnside's  method 


Also  contained  in  Figs.  4-4a,  4-5a,  and  4-6a  are  the  far-field 
patterns  of  the  monopole  as  if  it  were  on  an  infinite  ground  plane 
with  no  diffracted  field  contributions.  This  pattern  is  indicated 
by  the  + symbol  on  the  curve  and  the  pattern  with  the  edge 
diffraction  by  the  QD  symbol.  By  comparing  these  two  patterns,  the 
effect  of  the  edge  diffraction  is  quite  noticeable.  The  level 
differences  of  the  patterns  between  the  a and  b parts  of  Figs.  4-4, 
4-5,  and  4-6  is  because  the  pattern  including  diffraction  is 
normalized  with  the  pattern  from  the  monopole  on  an  infinite  ground 
plane.  Nonetheless,  the  similarity  between  the  patterns  calculated 
by  the  two  different  methods  is  apparent.  Also,  Fig.  4-6b  has  dis- 
continuities at  90°  and  270°,  which  are  similar  to  those  encountered 
with  the  hybrid  technique  before  the  second  order  edge  diffraction 
was  included. 


0.  Two  Monopoles  on  a Ground  Plane 

Computer-drawn  sketches  in  each  of  the  figures  in  this  section 
describe  the  geometry  for  that  particular  pattern.  An  explanation 
of  these  sketches  is  illustrated  with  an  example.  Fig.  4-7  is  a 
representation  of  the  computer  sketch  in  Fig.  4-8a.  Fourteen  inches 
is  the  corresponding  dimension  between  the  center  of  the  two  monopole 
array  and  the  edge  of  the  ground  plane,  while  0.2  inch  is  the  spacing 
between  the  tips  of  the  monopoles.  The  relative  size  of  the  monopoles 
in  the  computer  sketch  is  not  to  scale,  and  the  inclination  angle  of 
the  monopoles  to  the  ground  plane  indicated  by  the  sketch  does  not 
necessarily  represent  the  actual  angle  used.  The  angle  of  the 
monopoles  or  elements  with  respect  to  the  ground  plane  is 
indicated  to  the  right  of  the  computer  sketch  (20°,  in  this  case). 

All  patterns  in  this  section  begin  (0°)  on  the  underside  of  the 
ground  plane  so  that  endfire  for  the  array  occurs  at  180°  as 
indicated  in  Fig.  4-7  and  the  computer  sketches. 

Just  as  with  the  single  monopole  cases,  the  pattern  from  the 
antenna  without  the  edge  diffracted  fields  is  indicated  by  the  + 
symbol  in  the  figures  to  follow  and  the  total  field  with  diffraction 
by  the  QD  symbol.  By  comparing  these  two  patterns,  the  effects  of 
the  edge  diffraction  in  all  cases  can  be  studied. 

In  the  first  group  of  patterns.  Figs.  4-8,  4-9,  4-10,  three 
random  cases  are  chosen  and  compared  with  Burnside's  GTD  computer 
program  for  the  same  geometry. 

The  signal  levels  between  the  hybrid  method  and  Burnside's  GTD 
program  differ  because  the  field  including  diffraction  is  normalized 
with  the  field  containing  no  diffraction  for  comparison  purposes  in 
Figs.  4-3a,  4-9a,  4-10a.  The  discontinuities  occurring  in  some  of 
the  patterns  around  90  and  270  degrees  are  due  to  the  noninclusion 
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Fig.  4-7.  Representation  of  the  computer  drawing  below 
the  patterns  in  this  section. 


of  higher  order  diffraction.  It  is  not  worth  the  added  computational 
efforts  for  our  purposes  to  include  higher  order  diffraction  around 
the  edges  of  the  ground  plane.  The  two  programs  calculate  the 
patterns  by  different  methods  but  display  a great  deal  of  similarity. 

Figures  4-1  la,  b,  c are  the  patterns  of  the  monopoles  located 
symmetrically  on  the  ground  plane  for  800,  900,  and  1000  MHz.  The 
desired  pattern  shape  appears  to  be  obtainable  at  the  higher 
frequencies.  At  800  MHz,  the  null  at  180°  is  intolerable. 

Figures  4-1 2a,  b,  c and  Figs.  4-1 3a,  b,  c demonstrate  the  effects 
on  the  pattern  by  locating  the  monopoles  asymmetrically  on  the 
ground  plane.  The  monopole  or  element  spacing  and  angle  with  respect 
to  the  ground  plane  are  the  same  as  in  Figs.  4-1  la,  b,  c.  As  with 
the  symmetric  cases,  the  pattern  evolves  more  to  the  desired 
pattern  as  the  frequency  increases  to  1000  MHz.  As  the  frequency 
increases,  the  levels  of  the  nulls  in  the  region  of  main  concern 
from  120°  to  240°  grow  shallower. 

Figures  4-1 4a-f  display  the  effects  of  the  far-field  patterns 
due  to  changes  in  the  separation  between  the  two  monopoles  or 
elements.  As  the  spacing  is  increased,  the  pattern  shape  becomes 
less  desirable  due  to  the  deepening  of  the  center  null.  The 
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smaller  spacings  are,  therefore,  more  useful  in  achieving  the 
desired  pattern.  Thus,  0.3  inch  was  chosen  as  a reasonable 
spacing  for  the  remaining  patterns  in  the  section. 

Figures  4-15,  4-16,  and  4-17  illustrate  the  effects  on  the 
pattern  due  to  changes  in  the  angle  of  the  monopoles  with  respect 
to  the  ground  plane. 

Angles  less  than  30  degrees  were  excluded  because  the  close 
proximity  of  the  monopoles  to  the  ground  plane  in  the  true 
physical  situation  could  result  in  inadequate  input  impedances 
over  the  range  of  frequency.  Angles  greater  than  45  degrees  were 
also  tested  but  resulted  in  nondesirable  pattern  shapes. 

A monopole  angle  of  30  degrees  with  a tip  spacing  of  0.3  inch 
seems  to  yield  the  optimum  results.  In  Figs.  4-18  through  4-22 
these  geometrical  values  are  used  to  test  for  the  two  monopole  array 
at  different  locations  on  the  ground  plane. 

In  Figs.  4-18a,  b,  c the  discontinuities  around  70  to  110 
degrees  are  due  to  the  close  proximity  of  the  edge  to  the  monopole 
array.  Even  still,  the  general  pattern  shape  is  valid.  From 
studying  Figs.  4-18  through  4-22  a more  desirable  pattern  shape 
is  obtained  as  the  monopole  array  is  moved  closer  to  the  edge.  As 
the  array  is  shifted  toward  the  center,  deep  nulls,  particularly 
at  the  low  end  of  the  frequency  range  (800  MHz),  develop. 

Additional  patterns  are  contained  in  Appendices  I and  II. 
Appendix  I contains  patterns  with  the  center  of  the  monopole  array 
located  at  12  inches  from  one  ground  plane  edge  and  16  inches  from 
the  other  edge  with  the  monopoles  angles  set  at  30  and  45  degrees 
for  various  element  spacings. 

Appendix  II  has  patterns  with  the  center  of  the  monopole  array 
located  eight  inches  from  one  edge  and  20  inches  from  the  other. 

The  monopoles  angles  are  set  at  30  degrees,  and  different  element 
spacings  are  tested. 
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Fig.  4-8b.  Calculated  by  Burnside's  6TD  computer  program 
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Fia  4-10b . Calculated  by  Burnside's  computer  program. 
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Fig.  4-21 b. 


CHAPTER  V 

SUMMARY  AND  CONCLUSIONS 


This  research  and  development  effort,  which  investigates  the 
radiation  patterns  of  a proposed  UHF  satellite  antenna,  has 
utilized  a combination  of  two  powerful  computational  techniques, 
the  Method  of  Moments  and  the  Geometrical  Theory  of  Diffraction 
(GTD),  combined  into  a Hybrid  Technique.  The  Hybrid  Technique 
allows  us  to  investigate  the  performance  of  the  satellite  antenna 
at  various  locations  on  a 28"  by  28"  panel  without  making  an  assumption 
about  the  current  distribution  in  amplitude  or  phase  as  one  would 
normally  do  in  a pure  GTD  approach  to  the  problem.  Instead,  the 
current  distribution  is  treated  as  an  unknown  in  the  moment 
method  part  of  the  hybrid  technique. 

The  radiation  pattern  results  presented  in  Chapter  IV  have  been 
checked  against  an  independent  computer  program  and  against  measure- 
ments in  the  literature.  In  all  cases  agreement  is  good,  giving 
one  the  necessary  confidence  in  the  data  and  the  conclusions  that  can 
be  drawn  from  them.  It  is  surprisingly  apparent  from  the  data  that 
the  best  location  for  the  UHF  antenna  is  near  one  edge  of  the  panel , 
rather  than  near  the  center. 

There  are,  of  course,  other  antennas  and  scattering  objects 
on  the  panel  that  will  redirect  some  of  the  energy  and  cause  per-* 
turbations  in  the  radiation  patterns  from  those  shown  here.  The 
method  we  used  here  and  the  associated  computer  program  can  be 
extended  to  include  these  effects.  It  is  recommended  that  such 
an  effort  be  undertaken. 
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APPENDIX  I 


The  patterns  contained  in  this  appendix  are  those  of  a two 
monopole  array  with  the  array  center  located  12  inches  from  one 
ground  plane  edge  and  16  inches  from  the  opposite  edge. 

Figures  I-la-e  are  patterns  for  various  element  spacings 
at  800  MHz.  The  angle  of  the  elements  is  30  degrees.  Figures 
I-2a-e  and  Figs.  I-3a-e  have  the  same  geometric  set  ups  as  Figs. 
I-la-e  but  the  frequency  is  900  and  1000  MHz  respectively. 

Figures  I-4a-f,  I-5a-f,  and  I-6a-f  are  patterns  for 
various  element  spacings  with  the  angle  of  the  elements  set  at  45 
degrees.  Figures  I-4a-f  are  at  800  MHz  while  Figs.  I-5a-f  and 
Figs.  I-6a-f  are  at  900  and  1000  MHz  respectively. 
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Fig,  I-5e, 
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APPENDIX  II 


The  patterns  contained  in  this  appendix  are  those  of  a two 
monopole  array  with  the  array  center  located  at  8 inches  from 
one  ground  plane  edge  and  20  from  the  opposite  edge.  The  angle 
of  the  monopoles  with  respect  to  the  ground  plane  is  30  degrees. 
Figures  II-la-e  are  at  800  MHz  while  Figs.  II-2a-d  and  Figs, 
II-3a-e  are  at  900  and  1000  MHz  respectively. 
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